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ABSTRACT 

We use the Om statistic and the Genetic Algorithms (GA) in order to derive a nuU 
test on the spatiahy flat cosmological constant model ACDM. This is done in two 
steps: first, we apply the GA to the Constitution SNIa data in order to acquire a 
model independent reconstruction of the expansion history of the Universe H{z) and 
second, we use the reconstructed H{z) in conjunction with the Om statistic, which 
is constant only for the ACDM model, to derive our constraints. We find that while 
ACDM is consistent with the data at the 2cr level, some deviations from ACDM model 
at low redshifts can be accommodated. 

Key words: cosmology: theory - (cosmology:) dark energy - (cosmology:) cosmolog- 
ical parameters. 



1 INTRODUCTION 

In the previous decade it was discovered that the Universe 
is undergoing an accelerated expansion (Riess et al. 2004; 
Spergel et al. 2007; Readhead et al. 2004) . This acceleration 
is usually attributed either to a cosmic fluid with negative 
pressure dubbed Dark Energy or to an IR modification of 
gravity. In order to identify the properties of Dark Energy or 
the structure of the IR modification of gravity it is necessary 
to know to a high precision the rate of the expansion of the 
Universe, parameterized as H = ^ where a = is the 
scale factor and z is the redshift of the cosmological probe, 
as it measured by the observations. 

The behavior of the expansion of the Universe can be 
identified by studying two functions, the Equation of State 
(EoS) 'w{z) = — which can be rewritten as 



w{z) = -1 +-(1 



^ d\n{5H{zf) 
dz 



(1) 



where SH{zf = H{zf/Ho - Q.om{l + z)^ accounts for 
all terms in the Friedmann equation not related to matter 

and the deceleration parameter q{z) = ^ which can be 

rewritten as 



q{z) = -l+{l + z) 



dln{H{z)) 
dz 



(2) 



Obviously, the cosmological constant {w{z) — —1) corre- 
sponds to a constant dark energy density, while in general 
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w{z) can be time dependent. Also, an important param- 
eter is the value of the deceleration parameter today, ie 
q{z — 0) = go, which for the cosmological constant model in 
GR it is go = -H- 3noin/2. 

However, despite all the recent progress the origin of 
the accelerated expansion of the universe still remains un- 
known with many possibilities still remaining open, see for 
example Perivolaropoulos (2006). The simplest choice that 
agrees well with the data is a positive cosmological constant 
which has to be small enough to have started dominating the 
universe at late times. As it was demonstrated by the Seven- 
Year WMAP data (Komatsu et al. 2010), the cosmological 
constant remains the best candidate and has the advantage 
of having only one free parameter related to the properties 
of the Dark Energy. Nonetheless, this model fails to explain 
why the cosmological constant is so small that it can only 
dominate the universe at late times, a problem known as the 
coincidence problem and there are a few cosmological obser- 
vations which differ from its predictions (Perivolaropoulos 
2008; Perivolaropoulos & Shafieloo 2008). 

A very important complication in the investigation of 
the behavior of dark energy occurs due to the bias intro- 
duced by the parameterizations used. At the moment, there 
is a multitude of available phenomenological ansatze for the 
dark energy equation of state parameter w or dark energy 
density, each with its own merits and limitations (see Sahni 
& Starobinsky (2006) and references therein). The interpre- 
tation of the SNIa data has been shown to depend greatly 
on the type of parametrization used to perform a data fit 
(Sahni & Starobinsky 2006; Shafieloo, Sahni & Starobinsky 
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2009). Choosing a priori a model for dark energy can thus 
adversely affect the validity of the fitting method and lead 
to compromised or misleading results. 

The need to counteract this problem paved the way for 
the consideration of a complementary set of non-parametric 
reconstruction techniques (Daly & Djorgovski 2003; Wang 
& Mukherjee 2004; Saini 2003; Wang 2009; Clarkson & 
Zunckel) and model independent approaches (Sahni et al. 
2002; Alam et al. 2003; Shafieloo et al. 2006, Shafieloo 2007; 
Wang & Tegmark 2005; Shafieloo & Clarkson 2010). These 
try to minimize the ambiguity due to a possibly biased as- 
sumption for w by fitting the original datasets without using 
any parameters related to some specific model. The result 
of these methods can then be interpreted in the context of 
a dark energy model of choice. Non-parametric reconstruc- 
tions can thus corroborate parametric methods and provide 
more credibility. However, they too suffer from a different 
set of problems, mainly the need to resort to differentiation 
of noisy data, which can itself introduce great errors. 

In this paper, we present a method that can be used 
as a model independent approach in testing the standard 
cosmological model. This is done by using the Genetic Al- 
gorithms (GA) technique, first used in the analysis of SNIa 
data in Bogdanos & Nesseris (2009). The GAs represent a 
method for non-parametric reconstruction of the dark en- 
ergy equation of state parameter w, based on the notions 
of genetic algorithms and grammatical evolution. GAs are 
more useful and efficient than usual techniques when 

• The parameter space is very large, too complex or not 
enough understood, as is the case with dark energy. 

• Domain knowledge is scarce or expert knowledge is dif- 
ficult to encode to narrow the search space. 

• Traditional search methods give poor results or com- 
pletely fail. 

Naturally, therefore, they have been used with success in 
many fields where one of the above situations is encountered, 
like the computational science, engineering and economics. 
Recently, they have also been applied to study high energy 
physics (Becks, Hahn & Hemker 1994; AUanach, Grellscheid 
& Quevedo 2004; Rojo & Latorre 2004) gravitational wave 
detection (Crowder, Cornish & Reddinger 2006) and grav- 
itational lensing (Brewer & Lewis 2005). Since the nature 
of Dark Energy still remains a mystery, this makes it for us 
an ideal candidate to use the GAs as a means to analyze 
the SNIa data and extract model independent constraints 
on the behavior of the Dark Energy. In Section 2 we briefly 
describe the Om statistic while in Section 3 we provide an 
overview of the general methodology of the GA paradigm 
and flnally, we present our results in Section 4. 



2 THE OM STATISTIC 

The recently introduced Om diagnostic (Sahni, Shafleloo & 
Starobinsky 2008) (also look at Zunckel & Clarkson (2008)) 
enables us to distinguish ACDM from other dark energy 
models without directly involving the cosmic EOS. The Om 
diagnostic is deflned as: 




Figure 1. The Om diagnostic is shown as a function of red- 
shift for dark energy models with Qom = 0.27 and w = 
-1,-0.8,-0.6,-0.4,-0.2 (bottom to top). For Phantom mod- 
els (not shown) Om would have the opposite curvature. 



For dark energy with a constant equation of state w = const, 



h (x) = Q.QmX + (1 - f2om)a;° 



a = 3(1 + ii;) 



(4) 



(we assume that the universe is spatially flat for simplicity) . 
Consequently, 



Om{x) = Q,om + (1 - f^om)- 



(5) 



from where we find Om{x) = Q,om in ACDM, whereas 
Om{x) > Slom in quintessence {a > 0) while Om{x) < 
Qom in phantom (a < 0). We therefore conclude that: 
Om{x) — Q.om = if dark energy is a cosmological con- 
stant. Note that since SIa + f^Om — 1 in ACDM, this model 
contains SCDM (fiom = Ij^A = 0) as an important lim- 
iting case. Consequently the Om diagnostic cannot distin- 
guish between large and small values of the cosmological 
constant unless the value of the matter density is indepen- 
dently known. In other words, the Om diagnostic provides 
us with a null test of the cosmological constant. This is a 
simple consequence of the fact that h^ix) plotted against 
x'^ results in a straight line for ACDM, whose slope is given 
by Q,om- For other dark energy models the line describing 
Om{x) is curved, since the equality 



= constant 



(6) 



Om{x) 



h^ix) 



X — 1+z 



h{x) = H{x)/Ho .(3) 



(which always holds for ACDM for any x = 1 + z) is satis- 
fied in quintessence/phantom type models only at redshifts 
significantly greater than unity, when the effects of dark en- 
ergy on the expansion rate can safely be ignored. As a result 
the efficiency of the Om diagnostic improves at low redshifts 
(z < 2) precisely where there is likely to be an abundance 
of cosmological data in the coming years! 

For a constant EoS: l-|-ui ~ [Om{z) — flom]{l — ^0m)~^ 
at z <C 1, consequently a larger Om{z) is indicative of a 
larger w; while at high z, Om{z) Qom, as shown in Fig. [T] 

On the other hand, for quintessence as well as phantom 
the line describing Om(x) is curved, which helps distinguish 
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these models from ACDM even if the value of the matter 
density is not accurately known - see Fig. [T] 

In practice, the construction of Om requires a knowl- 
edge of the Hubble parameter, h{z), which can be deter- 
mined using a number of model independent approaches 
(Sahni et al. 2003; Alam et al. 2003; Shafieloo et al. 2006, 
Shafieloo 2007; Wang & Tegmark 2005; Shafieloo & Clark- 
son 2010). 

3 GENETIC ALGORITHMS 
3.1 Overview 

GAs were introduced as a computational analogy of adaptive 
systems. They are modelled loosely on the principles of the 
evolution via natural selection, employing a population of in- 
dividuals that undergo selection in the presence of variation- 
inducing operators such as mutation and crossover. The en- 
coding of the chromosomes is called the genome or genotype 
and is composed, in loose correspondence to actual DNA 
genomes, by a series of representative "genes". Depending 
on the problem, the genome can be a series (vector) of bi- 
nary numbers, decimal integers, machine precision integers 
or reals or more complex data structures. 

In order to evaluate the individual chromosomes a fit- 
ness function is used and reproductive success usually varies 
with fitness. The fitness function is a map between the gene 
sequence of the chromosomes (genotype) and a number of 
attributes (phenotype), directly related to the properties of 
a wanted solution. Often the fitness function is used to de- 
termine a "distance" of a candidate solution from the true 
one. Various distance measures can be used for this purpose 
(Euclidean distance, Manhattan, Mahalanobis etc.). 

The algorithm begins with an initial population of can- 
didate solutions, which is usually randomly generated. Al- 
though GA's are relatively insensitive to initial conditions, 
i.e. the population we start from is not very significant, but 
using some prescription for producing this seed generation 
can affect the speed of convergence. In each successive step, 
the fitness functions for the chromosomes of the popula- 
tion are evaluated and a number of genetic operators (mu- 
tation and crossover) are applied to produce the next gener- 
ation. This process continues until some termination criteria 
is reached, e.g. obtain a solution with fitness greater than 
some predefined threshold or reach a maximum number of 
generations. The later is imposed as a condition to ensure 
that the algorithm terminates even if we cannot get the de- 
sired level of fitness. 

The various steps of the algorithm can be summarized 
as follows: 

(i) Randomly generate an initial population A/(0) 

(ii) Compute and save the fitness for each individual m 
in the current population M{t). 

(iii) Define selection probabilities p(m) for each individ- 
ual m in M(t) so that p{m) is proportional to the fitness. 

(iv) Generate M{t-\-l) by probabilistically selecting indi- 
viduals from M(t) to produce offspring via genetic operators 
(crossover and mutation). 

(v) Repeat step 2 until a satisfying solution is obtained, 
or a maximum number of generations reached. 



only depend on the choice of the grammar and therefore it 
can only affect on how fast the algorithm will converge to 
the minimum. Using the wrong grammar may result in the 
GA being trapped in a local minimum. Also, two important 
parameters that affect the transition from the population 
M{t) to M{t -\- 1) are the selection rate and the mutation 
rate. The selection rate is typically of the order of 10% and 
it affects how many of the individuals, after they are ranked 
with respect to their fitness, will be allowed to produce off- 
spring. The mutation rate is usually of the order of 5% and 
it expresses the probability that an arbitrary part of the ge- 
netic sequence will be changed from its previous state. This 
is done in order to maintain the genetic diversity from one 
generation of a population to the next. 

The paradigm of GAs described above is usually the one 
applied to solving most of the problems presented to GAs. 
Though it might not find the best solution, more often than 
not, it would come up with a partially optimal solution. A 
more detailed overview and its application in cosmology can 
be found in Bogdanos & Nesseris (2009). 

The difference of the GA over the the standard analysis 
of the data, ie having an a priori defined model with some 
free parameters, is that the obtained values of the best-fit 
parameters are model-dependent and in general models with 
more parameters tend to give better fits to the data. This is 
where the GA approach starts to depart from the ordinary 
parametric method. Our goal is to minimize a function, not 
using a candidate model function for the distance modulus 
and varying parameters, but through a stochastic process 
based on a GA evolution. This way, no prior knowledge of 
a dark energy model is needed to obtain a solution and our 
result will be completely parameter-free. 

In simpler terms, the GA method does not require an a 
priori assumption for a DE model, but uses the data them- 
selves to find this model. Also, it is parameter free as the end 
result does not have any free parameters that can be changed 
in order to fit the data. So, in this sense this method has far 
less bias any of the standard methods for the reconstruction 
of the expansion history of the Universe. This is the main 
reason for the use of the GAs in this paper. 



3.2 General Methodology 

We first outline the course of action we follow to apply the 
GA paradigm in the case of SNIa data. For the application 
of GA and grammatical evolution (GE) on the dataset, we 
use a modified version of the GDF (Tsoulos et al. 2006) tool 
0, which uses GE as a method to fit datasets of arbitrary 
size and dimensionality. This program uses the tournament 
selection method for crossover. GDF requires a set of input 
data (train set), an (optional) test sample and a grammar 
to be used for the generation of functional expressions. The 
output is an expression for the function which best fits the 
train set data. 

Since the SNIa datapoints are given in terms of the 
distance modulus fiobs{zi), the fitness function for the GA 
we have chosen is equal to — xliv/m where 



We should stress that the initial population M(0) will 



^ |http://cpc.cs. qub .ac.uk/ summaries / ADXC] 
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2 _ {fJ-obs{Zi) - llGAiZr)f , - 

XSNIa — 2 ' \') 

z — 1 ^ 

where for the Constitution set = 397 and ^ are the er- 
rors due to flux uncertainties, intrinsic dispersion of SNIa 
absolute magnitude and peculiar velocity dispersion and 
liGA{zi) is the reduced distance modulus obtained for each 
chromosome of the population by the GA. 

An advantage of our method is the fact that neither the 
xliv/a nor the ^ga{z) depend on any parameters. Also, as 
it can be seen from eq. ((7|, our method does not make an 
assumption of flatness or not in order to fit the data. Flat- 
ness (or non-flatness) only comes into play when one tries to 
find the luminosity distance dL(z) and the underlying dark 
energy model given the best-fit form of ^ga{z). So, in this 
aspect the result produced by the GA is independent of the 
assumption of flatness as well. 

The GA evaluates (O in each evolutionary step for every 
chromosome of the population. The one with the best fitness 
will consequently have the smallest x|jv/a s^nd will be the 
best candidate solution in its generation. Of all the steps in 
the execution of the algorithm, the evaluation of the fitness 
is the most expensive. GDF is appropriately modified to use 
((TJ as the basis of its fitness calculation. 

After the execution of the GA, we obtain an expression 
IJ-ga{z) for the reduced distance modulus as the solution of 
best fitness and the corresponding XsNia- Using this we can, 
through differentiation, obtain other functions or parame- 
ters of interest, such as the Om statistic or the deceleration 
parameter q{z). For example, for a flat Universe the Hubble 
parameter will be given by 

and then the deceleration parameter can be found by Eq. ([2]). 

A possible complication of this method is that it gives 
no direct way to estimate the errors for the derived param- 
eters. One cannot expect to get an estimate by just running 
the algorithm many times and obtaining slightly different 
parameters. The GA usually tends to converge at the same 
solution for a given dataset, unless we change significantly 
the population size or the number of generations. A way to 
circumvent this problem is to use a bootstrap Monte Carlo 
simulation to produce synthetic datasets and rerun the al- 
gorithm on them. We can thus obtain a statistical sample of 
parameter values which will allow us to estimate the error. 

We sketch the procedure we will follow Press et al. 
(1992): 

(i) The GA is applied on the original SNla dataset with 
the chosen execution parameters and a solution for jlGA(z) 
is obtained. 

(ii) Generate a number of synthetic datasets by drawing 
each time the same number of data points with replacement 
from the original set. 

(iii) The GA is rerun for the synthetic datasets. 

(iv) A new set of values for the Om statistic and the de- 
celeration parameter q{z) is generated. 

(v) The 95% error of the desired parameter can be found 
by taking the 2.5 and the 97.5 percentiles of the bootstrap 
distribution (Efron 1982). 



Using the above steps, we can obtain error estimates for any 
desired parameter. 

3.3 An example 

In this Section we will briefly describe a simple exampl43 of 
how he GA determines the best-flt. We will avoid describ- 
ing the technicalities, like the binary representation of the 
solutions, and instead we will concentrate on how the gen- 
erations evolve, how many and which individuals are chosen 
for the next generation etc. 

As we mentioned earlier, the choice of the grammar is 
very important, however in order to keep our example as 
simple as possible we will assume that our grammar includes 
only basic functions like polynomials x, etc, the trigono- 
metric functions sin{x), cos{x), the exponential and the 
logarithm ln{x). 

The first step in the GA is setting up a random initial 
population A/(0) which can be any simple combination of 
these functions, eg fj,GA,i{z) — ln{z), fJ,GA,2{z) = —1 + z + z^ 
and i-iGA,3{z) = sin{z). The number of number of candidate 
solutions (chromosomes) in the genetic population is usually 
a few hundreds and later on we will use the value of 500. 

Next, the algorithm measures the fitness of each so- 
lution by calculating their for our simple example 
xl = 80579.9, xi = 292767.0 and Xs = 412928.0. The se- 
lection per se is done by implementing the "Tournament 
selection" method which involves running several "tourna- 
ments", sorting the population with respect to the fitness 
each individual and after that a fixed percentage, adjusted 
by the selection rate (see the previous subsection), of the 
population is chosen. As we mentioned, the selection rate is 
of the order of 10% of the total population, but lets assume 
that the two out of the three candidate solutions (^ga,i(z) 
and jJ.GA,2{z)) are chosen for the sake of simplicity. 

The reproduction of these two solutions will be done by 
the crossover and mutation operations. The crossover will 
randomly combine parts of the "parent" solutions, for exam- 
ple in one such realization this may be schematically shown 
as 

fJ-GA,l{z) ® flGAaiz) {fJ-GA,l{z),flGA,2{z),fJ-GA,:i{z)) 

After this is done, the GA will proceed to implement 
the mutation operation. The probability for mutation is ad- 
justed by the mutation rate, which as we mentioned is typi- 
cally of the order of 5%. In our example this may be a change 
in the power of some term or the change in the value of a 
coefficient. For example, for the candidate solution flGA,3{z) 
this can be schematically shown as fiGA,3{z) = —l-\-ln{z) — >■ 
— l+ln{z'^), where the power of the z term was mutated from 
1 to 3. 

Finally, at the end of the first round we have the three 
candidate solutions 

M(l) = {fjLGA,l(z),flGA,2{z), flGA/iiz)) = 

{ln{z^), -1 + ln{z^), -1 + ln{z^)) 

^ This is only a schematic description of how the GA works and 
this is done solely for the sake of explaining the basic mechanisms 
of the GA. 
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At the beginning of the next round the fitness of each indi- 
vidual will again be determined, which for our population 
will be (xi,X2,X3) = (3909.2,18435,113327), and the se- 
lection and the other operations will proceed as before. It 
is easily seen that even after one generation the oi the 
candidate solutions has decreased dramatically and as we 
will see in the next section, it only takes a bit more than 
100 generations to reach an acceptable a-nd about 500 to 
1000 generations to start converging on the minimum. 

After a predetermined number of generations, usually 
on the order of 10^, has been reached or some other criteria 
have been met, the GA will finish, with the best candidate 
solution of the last generation being the best-fit to the data. 



Table 1. The different cases considered in the analysis. The Ba- 
sic Grammar includes polynomials, trigonometric functions, ex- 
ponentials and logarithms. For ACDM the best-fit corresponds to 
Um = 0.289. 



Case 1 


Basic grammar 


\^ ■ 


= 465.19 


Case 2 


Basic grammar -|- Legendre polyn. 


y^ ■ 


= 477.87 


Case 3 


Only polyn. in grammar (3000 gen.) 


y^ ■ 


= 468.19 


Case 4 


Only polyn. in grammar (6000 gen.) 


y^ ■ 


= 462.56 


ACDM 


H{z)'^ = H2(Q„(1 + z)3 + i^n^) 


y^ - 


= 465.51 



4 RESULTS 

We will now examine a number of different configurations, 
highlighted in Table [T] In order to keep our analysis as 
simple as possible we will use a Basic Grammar that in- 
cludes polynomials like x, etc, the trigonometric functions 
sm{x), cos{x), the exponential and the logarithm ln{x) 
and a second grammar that includes the Basic Grammar 
and the Legendre polynomials. For the analysis we will use 
the Constitution SNIa dataset of Hicken et al. (2009) con- 
sisting of 397 SNIa. The steps followed for the usual min- 
imization of ((7| in terms of its parameters are described 
in detail in Nesseris & Perivolaropoulos (2004); Nesseris & 
Perivolaropoulos (2005); Lazkoz, Nesseris & Perivolaropou- 
los (2005), Nesseris & Perivolaropoulos (2007). 

In Fig. [2] we show the fitness (equal to — x^) as a func- 
tion of the number of generations. Clearly, all configurations 
seem to have to have converged very early or to be very 
close to converging to their minimum. An exception is Case 
3 which even after 3000 generations still has not converged 
and we had to run the simulation for another 3000 gener- 
ations to check its convergence. Since we want to compare 
the configurations with equal standards, ie number of gen- 
erations etc, we will treat this case (Case 4) separately. 

As it can be seen from Table [1] and Fig. [3] most config- 
urations agree well with the data and actually have a mini- 
mum x^ comparable or even better with respect to ACDM. 
An exception to this is Case 2, where the best fit differs from 
ACDM by a Ay^ — 12 thus providing a much poorer fit. The 
reason for this is the fact that by increasing the grammar 
we effectively increase the landscape in which the algorithm 
has to search and therefore, we are making the fit poorer as 
the GA will take much much longer to converge. 

In Figs. [4] and [5] we show the corresponding results 
of each case for the Om statistic and the deceleration pa- 
rameter q{z) (our best-fit function corresponds to the cyan 
line). Due to the presence of trigonometric functions into 
the best fit and the fact that the Om statistic is derived 
from the Hubble parameter by differentiation, Case 1 ex- 
hibits strong oscillations. This is also apparent in the cor- 
responding plot for the deceleration parameter for Case 2 
(magenta line. Fig. [S]). 

As it can be seen in Fig. [T] a constant w (but differ- 
ent from w = —1) forces the Om statistic to be a monotonic 
function that is located completely either above or below the 
constant case (Om — 0.289). However, none of the curves 
in Fig. |4] are monotonic functions and all of them cross the 




Figure 2. The fitness (equal to — y^) eis a function of the number 
of generations. The magenta line corresponds to case 1, the green 
line to case 2 and the blue line to case 3. The inset graphic shows 
case 4. 



constant Om case several times. Also, all cases remain rea- 
sonably close to the ACDM value of Sim = 0.289 for in- 
termediate redshifts but deviate from that value at high or 
low redshifts. Therefore, this could be the effect of a time- 
dependent equation of state w{z) instead of a constant w. 
Comparing with Fig. [1] it is possible that this behaviour is 
due to an equation of state w{z) more negative than -1. 




Figure 3. The SNIa distance modulus against the redshift. The 
red line represents the best-fit ACDM model with Qm = 0.289 
while the magenta, blue, green and cyan (best-fit) lines corre- 
spond to cases 1-4 respectively. 
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This is also supported from Fig. |4] and Fig. [5] as at low 
redshifts all cases predict a much different forms of accel- 
eration than the best fit ACDM. In Fig. 0] and Fig. [5] our 
best fit results show an overall slowing down of the accel- 
eration at the low redshifts (which is in concordance with 
results presented in Shafieloo et al. (2009) with a very rapid 
increase of acceleration at the very low redshifts. These are 
very similar with the results presented in Fig. 2 of Shafieloo 
& Clarkson (2010) where smoothing method of reconstruc- 
tion of the expansion history of the universe is being used. 
It is obvious that when two very different model indepen- 
dent approaches reach to a similar result, we must be on 
the right track and very close to the best possible result. 
However, theoretical interpretation of these results is be- 
yond the scope of this paper. In fact there might be very 
different mechanisms that result to a similar expansion his- 
tory of the universe as we reconstructed in this paper (e.g 
see Csaki, Kaloper, and Terning (2002); Alam et al. (2003), 
Kelly et al. (2010); Sullivan et al. (2010)). 

In Fig. [6] we show the Om statistic as a function of the 
redshift z. The black line corresponds to the best fit of case 
4, while the gray-shaded area to the 2a error region. The 
error region was calculated by implementing a bootstrap 
monte-carlo simulation as discussed in the previous section. 
Note that the overall slowing down of the acceleration at 
the low redshifts cannot be clearly seen from this figure as 
the errors are quite symmetrical around the best fit ACDM 
value. However, as we mentioned earlier this is clearly seen 
by the behavior of the deceleration parameter in Fig. [5] 

As it can be seen in Fig. (6] ACDM remains consistent 
with the data at the 2a level but at the same time, some 
deviations from the standard ACDM model can be accom- 
modated. Due to the large errorbars of the current datasets, 
many dark energy models are still consistent with them but 
we expect this to change in the near future with the advent 
of high quality SnIa observations. 

In Fig. [7] we show a histogram of the bootstrap dis- 
tribution found from the monte-carlo simulation that was 
used to create the error regions. Clearly, the distribution is 
quite symmetrical and centered, so we are confident that our 
error region in Fig. [6] was correctly reconstructed. Finally, 
we also tested the GA on a mock ACDM dataset and we 
found that it can reconstruct the luminosity distance curve 
to within 0.1%, so we are confident that is does not overfit 
the data, but instead it can successfully detect the underly- 
ing DE model. 



5 CONCLUSIONS 

We used the Om statistic and the GAs in order to derive 
a null test on the cosmological constant model ACDM. Our 
interest in the GAs stems from the fact that they represent 
a method for non-parametric reconstruction of the expan- 
sion history of the universe, based on the notions of genetic 
algorithms and grammatical evolution. These kinds of algo- 
rithms are more useful and efficient than usual techniques 
especially when the problem under examination is not well 
understood, as is the case with dark energy. Since the nature 
of dark energy still remains a mystery, this makes it for us 
an ideal candidate to use the GAs as a means to analyze 
the SNIa data and extract model independent constraints 
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Figure 4. The Om statistic against the redshift. The black 
dashed line represents the ACDM value of Qm = 0.289 while 
the magenta, blue, green and cyan (best-fit) lines correspond to 
cases 1-4 respectively. 
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Figure 5. The deceleration parameter q against the redshift. The 

red line represents the best-fit ACDM model with Vlrn = 0.289 
while the magenta, blue, green and cyan (best-fit) lines corre- 
spond to cases 1-4 respectively. 
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Figure 6. The Om statistic as a function of redshift. The black 
line corresponds to case 4 (the best-fit), while the gray-shaded 
area to the 2a error region. The error region was calculated by 
implementing a bootstrap monte-carlo simulation (see text for 
details). 
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Figure 7. Histogram of the bootstrap distribution found from the 
monte-carlo simulation that was used to create the error regions 
of Fig.E] 



on the behavior of the Dark Energy. On the other hand, the 
Om diagnostic (Sahni et al. 2008) enables us to distinguish 
ACDM from other dark energy models without directly in- 
volving the cosmic EoS. 

Our methodology is completely model independent and 
it can be summarized in two steps: first, we applied the GA 
to the Constitution SNIa data in order to acquire a model 
independent reconstruction of the expansion history of the 
Universe H(z). After we reconstructed H{z), the second step 
was to use it in conjunction with the Om statistic and derive 
our null test. 

Our main results was that ACDM remains consistent 
with the data at the 2a level, see Fig. [6] but at the same 
time, some deviations from the standard ACDM model can 
be accommodated and this is also in accordance with the 
data especially at low redshifts (see the behavior of the best 
fit in Fig. El). 

However, we should mention that the slowing down at 
low redshifts for the Constitution set mentioned earlier can 
also be seen with the standard analysis with the CPL ansatz 
(see Sanchez, Nesseris & Perivolaropoulos (2009)) but only 
for the Constitution set. For the other SNIa datasets the 
reverse trend (speeding up the acceleration at low z) was 
observed. For this reason and due to the fact that the current 
data still have quite large errors we believe that the full 
potential of our method, while it is very promising, will only 
be realized in the near future when more high quality SNIa 
data will become available. When this happens, the error 
region at low z in Fig. [6] will be small enough for the Om 
statistic to completely discriminate between ACDM and the 
various Dark Energy models. 
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